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ABSTRACT 

Context. There are many unknowns in the formation of subdwarf B stars. Different formation channels are considered to be possible 
and to lead to a variety of helium-burning subdwarfs. All seismic models to date, however, assume that a subdwarf B star is a post- 
helium-flash-core surrounded by a thin inert layer of hydrogen. 

Aims. We examine an alternative formation channel, in which the subdwarf B star originates from a massive (> ~2 M Q ) red giant 
with a non-degenerate helium-core. Although these subdwarfs may evolve through the same region of the log g - T sff diagram as 
the canonical post-flash subdwarfs, their interior structure is rather different. We examine how this difference affects their pulsation 
modes and whether it can be observed. 

Methods. Using detailed stellar evolution calculations we construct subdwarf B models from both formation channels. The iron 
accumulation in the driving region due to diffusion, which causes the excitation of the modes, is approximated by a Gaussian function. 
The pulsation modes and frequencies are calculated with a non-adiabatic pulsation code. 

Results. A detailed comparison of two subdwarf B models from different channels, but with the same logg and r cff , shows that their 
mode excitation is different. The excited frequencies are lower for the post-flash than for the post-non-degenerate subdwarf B star. 
This is mainly due to the differing chemical composition of the stellar envelope. A more general comparison between two grids of 
models shows that the excited frequencies of most post-non-degenerate subdwarfs cannot be well-matched with the frequencies of 
post-flash subdwarfs. In the rare event that an acceptable seismic match is found, additional information, such as mode identification 
and logg and T e g determinations, allows us to distinguish between the two formation channels. 

Key words, subdwarfs - stars: evolution - stars: oscillation - methods: numerical 



1. Introduction 

Commonly, subdwarf B (sdB) stars are identified as extreme hor- 
izontal branch (EHB) stars, and they are believed to be post-He- 
core-flash products with core masses ~0.5 M Q surrounded by 
a very thin inert H-envelope ( |Heber|[T986} |Saffer et aL"1|1994[ l. 
From a single stellar evolution point of view, this can be ex- 
plained by enhanced mass loss of stars close to He-ignition with 
very lightly bound envelopes ( jD'Cruz et al.|1996 1, i.e. stars with 
degenerate cores near the tip of the red giant branch (RGB). 
However, as they are frequently observed in binaries (e.g. Allard 
|et al.[ 1994; Morales -Rueda et al.|20 06 ), binary interactions most 
likely play an important role in their formation. |Han et aL ( 2002 ) 
explored the main binary evolution channels that can produce 
sdB stars: common-envelope ejection (CEE), stable Roche lobe 
overflow (RLOF), and helium white dwarf mergers. They found 
that the sdB mass distribution may be much broader than previ- 
ously thought, 0.3 - 0.8 M instead of 0.4 - 0.5 M Q . The sdB 
stars with non-canonical masses follow from mergers or massive 
(> ~2 M Q ) progenitors that ignite helium quiescently, where the 
latter can be a subchannel of either the CEE channel or the sta- 
ble RLOF channel. Binary population synthesis shows that the 
massive progenitors do not contribute significantly to the sdB 
population ( |Han et al.|2003) >. But one should keep in mind that it 



is asumed in such studies that CE evolution is described by the 
a formalism, i.e. that the CE ejection is driven by the orbital en- 
ergy. Because the physics of the CE phase is poorly understood, 
other scenarios should not be excluded a priori. For example, 
the y-formalism proposed by Nelema ns et al.| ([2000 ), based on 
the angular momentum equation rather than the energy equation, 
provides an alternative description. In this case, the massive red 
giants cannot be ruled out as possible progenitors of post-CE sdB 
stars ( Hu et al.|2007 1. We therefore want to explore the possibil- 
ity of this neglected class of progenitors in a different manner, 
by using the seismic properties that have been observed in some 
sdB stars. 

Although the post-flash and the post-non-degenerate sdB 
stars can appear in the same logg - T e s region, their interior 
structure is quite different. In particular, the chemical composi- 
tion profiles differ greatly depending on whether helium ignited 
in a flash or quiescently. For example, the canonical post-He- 
flash sdB star has a very narrow He-H transition zone, while 
the sdB star created from a more massive progentior has a much 
broader H-profile. This is a direct result of the differing chemi- 
cal compositions between low-mass and high-mass stars on the 
RGB, owing to the different convective regions during the main- 
sequence and RGB evolution. We examine whether this differ- 
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ence in the interior structure will result in observable differences 
in the pulsation modes. 

The sdB pulsators consist of two classes, the short-period 
variable EC 14026 stars ( |Kilkenny et al. 1997), and the long- 
period variable PG 1716 stars ( |Green et al. 2003 ). The rapid os- 
cillations in EC 14026 stars are interpreted in terms of low-order 



p-modes (Charpinet et al. 19961, driven by the /e-mechanism op- 
erating in the iron opacity bump. The same mechanism has been 
shown to excite long-period, high-order g-modes in the cooler 
models (Fontai ne et al.|20 03). The local iron enhancement nec- 
essary in the driving region around log T « 5.3 is due to the 
competing diffusion processes of radiative levitation and gravita- 
tional settling. It is well-known that the opacities play an impor- 
tant role in the study of the pulsations. Seaton & Badnell (2004) 
showed that the iron opacity bump is situated at slightly higher 
temperatures using OP opacities ( Seat on et al.|1994[ [Badnell & 
Seaton|2003[) compared with OPAL opacities (Iglesias & Rogers 



1996). ]Jeffery & Sai o (2006) found that, using OP opacities and 
nickel enhancement in addition to iron, the theoretical instabil- 
ity strip of g-mode sdB oscillators is more consistent with obser- 
vations. For our purposes it is sufficient to use OPAL opacities 
and iron enhancement, since we are interested in the relative dif- 
ferences between two types of sdB stars. We acknowledge the 
importance of including the effect of OP opacities and nickel 
enhancement in further detailed studies. 

The details of the computations are given in ^2] The results 
are presented in ^3] In 4 3.3 we compare the detailed physical 
characteristics of two reference models with different formation 
histories. In § |3.4| we compare the frequency characteristics glob- 
ally between two grids of models. The results and conclusions 
are discussed in $4l 



2. Computations 

2.1. The evolution calculations 

We constructed sdB structure models with the stellar evolution 
code developed by Eggleton ( 197l]|1972||1973[), Eggleton et al.| 



1973] ), and updatedby THan et al.| ( fT994| > and |Pols et al.| ( |1995 



The updated version of the code uses an equation of state 
that includes pressure ionization and Coulomb interaction, nu- 
clear reaction rates from Caughla n et al.| ( 1985 i and Caughlan 
& Fowler| ( |l988| ), and neutrino loss rates from |Itoh et al.| (|1989 
1992 ). Both convective and semi-convective mixing are treated 
as diffusion processes. It is assumed that mixing occurs in re- 
gions where 



V rad > V ad - <W(2.5 + 20/3 + 16/3 1 ), 



(1) 



where f3 is the ratio of radiation pressure to gas pressure and 6 ov 
is the overshooting parameter. Schroder et al. ( 1997 1 showed that 
6 m =0.12 gives the best fit to observations of £ Aurigae binaries, 
which corresponds to an overshooting length of ~0.25H p . For 
our comparative study, it suffices to adopt <S 0V = 0.12, but keep 
in mind that core overshooting can in fact also be probed by 
asteroseismology, e.g. |Aerts et al. ( 2003| ). 

We evolved stars assuming a chemical composition of X — 
0.70 and Z = 0.02. We used a mixing-length parameter (the 
ratio of the mixing-length to the local pressure scaleheight) of 
a — I /H p — 2.0. If not mentioned otherwise, the opacity tables 



were constructed by combining the OPAL opacities (Iglesias & 
Roge rs|1996[) with t he conductive opacities (Hubbard & Lampe 
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code bylEldridge & Tout 



1975[), as implemented in the Eggleton 
p004) . 



We started by evolving zero-age main-sequence (ZAMS) 
models in the range 1-3 Mq to the tip of the RGB, adopting 
a Reimer's mass loss rate (R eimers|1975| ), 



= 4x 10 _13 77 



(R/R G )(L/L Q ) 
(M/Mq) 



[M yr 



(2) 



with an efficiency of rj — 0.4 ( Iben & Renzini 1983 Carraro et al. 
1996). For simplicity, we did not include mass-loss on the EHB. 
|Unglaub & Bues| ( 2001| ) showed that, if the observed chemical 
abundances are the result of the combined effects of diffusion 
and mass loss, the sdB mass-loss rate should be in the range 
10~ 14 < M(M yr~') < 10~ 12 . T his is consistent with the rates 
found by |Vink & Ca ssisi (2002 ) for radiation-driven wind mod- 
els. They also showed that these rates are too low to have a direct 
effect on the sdB evolution. 

At the RGB tip, we removed the envelope, while keeping 
the chemical compositions fixed. Thus, we assume that the mass 
transfer happens on a much shorter timescale than the nuclear 
timescale. This is a reasonable assumption for sdB stars in short- 
period binaries formed by CE ejection, which is the majority of 
the observed sdB stars (Max ted et al.|200T| ) and the focus of our 
study here. In the case that the He-flash occurs, zero-age hor- 
izontal branch models were artificialy created from a 2.25 M Q 
He-core-burning star, where we reset the chemical compositions 
to the values before the flash. This treatment is not rigorously 
valid. Full evolutionary models of the He-flash show that the C 
abundance in the He-core can increase up to ~5% (Piersanti et al. 
2004 Serenelli & Weiss|2005[ ). Since the p-modes are not sen- 
sitive to the core, we are not worried about this. An interesting 
scenario is an sdB star that is formed by a late He-core flash on 
the white dwarf cooling curve ( |CasteUani~& Castellani 1993}. 
In such a case the He-flash-driven convection zone can pene- 
trate into the H-rich layers, resulting in a surface enrichment of 
He an d C ( |Brown et al.|200lHSchTattl et al.|200"T]|Cassisi et al.| 
2003). We note that this might influence the pulsations, but we 



will not discuss this scenario further here. 

On the EHB, we used for temperatures log T > 7 the 
same opacities as mentioned above. In the outer layers of the 
star, log T < 7, where the pulsation driving region is located, 
the opacities were calculated by interpolating between several 
OPAL tables computed with iron abundance enhanced by fac- 
tors of / = 1, 2, 5 , and 10 relative to solar, thus X(Fe) = 
0.071794Z/. The abundances of the other heavy elements are 
decreased such that the overall metallicity is kept constant as in 
IMiglio et al.| < |2007) . 

2.2. The oscillation calculations and iron accumulation 

We adapted the Eggleton evolution code so that the output is 
suitable for pulsation calculations. In practice, this implied cal- 
culating some additional physical quantities during the evolu- 
tion, and modifying the mesh to have sufficient meshpoints in 
the stellar envelope. The seismic properties of the stellar mod- 
els are then calculated with two pulsation codes. The adiabatic 
code OSC by Scufl aire et al.| ( [2007] ) is used to obtain the approx- 
imate frequencies, which are used as a first guess in the linear 
non-adiabatic code MAD by Dupret (2001 ). We determined the 
theoretical frequency spectrum up to / = 2, since it is expected 
that higher order modes are geometrically cancelled. Charpinet 
|et al.| ( fT996| ) has established that the excitation of sdB oscilla- 
tons is related to a local enrichment of iron in the stellar enve- 
lope caused by diffusion. Radiative levitation is expected to set 
up significant chemical gradients within a diffusion timescale of 
~10 5 yr, and consequently iron accumulates around log T ~ 5.3 
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(Michaud et al. 1989 1 Chayer et al. 1995). Time-dependent diffu- 
sion calculations |Fontaine et al.|20 06 ) show that, after ~10 5 yr, 
many pulsation modes are excited. Since element diffusion is not 
treated in the evolution code, we used an approximation for the 
iron accumulation, assuming that the iron only affects the stellar 
structure through the opacity. At each timestep of the evolution 
calculations, the iron enhancement factor / is increased with a 
Gaussian centered at log T = 5.3, 



df 
dt 



(W/10) 3 



exp( 



(log T- 5.3Y 



(3) 



with the initial condition f(t = 0) — 1. The width <x 2 = 0.05 
and accumulation timescale t = 4 x 10 5 yr are chosen such 
that iron is only increased in the region 4.5 < log T < 6.1, and 
lim^co/ = 10, which is loosely based on the time-dependent 
diffusion calculations of Fontaine et al.| ( |2006| l and the equilib- 
rium profiles of Charpin et et al.| ( |1997) . Our parametric approxi- 
mation is rather ad hoc, but since we are interested in the relative 
differences between two different scenarios, the exact shape of 
the iron profile is not crucial here. We will discuss the influence 
of the iron profile on the pulsations and the evolution of the star 
in CTT] 



3. Results 

3.1. Effects of the iron accumulation 

In Fig.[T] we show / throughout the star for different ages of an 
sdB star. Note that the temperature range 4.5 < log T < 6.1 cor- 
responds to a very narrow mass shell of ~10~ 6 M Q . In Fig. [I]:, we 
included the backreaction of convective mixing on our paramet- 
ric iron profile during the evolution. Note that the iron abundance 
is homogeneous near log T = 5.3 (log q = -10) and log T = 4.6 
(log q - -12.5). This is caused by two narrow convective layers 
due to iron and helium ionization, respectively. Interestingly, we 
found that the convective region around log T — 5.3 would not 
be present without iron accumulation. We determined that the 
slightly perturbed iron profile has a negligible effect on the driv- 
ing and the pulsation frequencies. Moreover, since our descrip- 
tion of iron accumulation is approximate, we have not included 
this effect in subsequent calculations. 

In Fig. |2j we show the effect of different iron abundance 
profiles on the sdB evolution in the logg - T e ff diagram. We 
compared our parametric approach (Eq. [3]) with the case of no 
iron enhancement, and a uniform enhancement of / = 10 in the 
whole envelope, as used in studies of mode excitation (Jef fery &| 
|Saio||2006| . It is evident that using / = 10 influences the evo- 
lution drastically, while our Gaussian parametrization of / has 
little effect. Higher iron abundances indeed give higher opaci- 
ties, and thus larger stellar radii. With Eq. [3] we only increase 
iron in a relatively small region, resulting in a minimal effect on 
the stellar structure. However, a minimal change of the stellar 
structure can result in a visible shift of the pulsation frequencies, 
as shown by Fontaine et al. (2006). 

We also compared the effect on the excitation, and find that 
our approximation of the iron profile can excite almost as many 
modes as the / = 10 enhancement, see Fig. [3^-c This is un- 
derstood in terms of the driving mechanism in sdB stars, which 
is associated with the iron opacity bump such that accumulat- 
ing iron in this driving region is sufficient for the excitation of 
the pulsation modes, see Fig. |3jl-f. Thus, with our parametric 
approach to iron accumulation, the issue of excitation can be ad- 
dressed while keeping the effect of the iron profile on the stellar 




Fig. 1. The iron enhancement factor / throughout the star as a 
function of (a) the outer mass fraction: \ogq = log(l - jj^J, 
and (b) the temperature. In panel (c), we show the effect of con- 
vective mixing on our parametric iron profile. The sdB star has 
M COK = 0.47 M and M env o = 10" 4 M , and was constructed 
from a 1.00 M ZAMS model. The profiles from bottom to top 
correspond to sdB ages 10 5 , 10 6 , 10 7 , 1.8 x 10 8 yr, where the last 
model is at the end of core-He-burning. 



structure realistic. Furthermore, these are the first evolutionary 
models of sdB stars that include the effect of iron accumulation, 
albeit in an approximative manner. 

3.2. The stellar models 

Following the procedure as described in ^2] we constructed 
a grid of canonical, i.e. post-He-flash, sdB models (hereafter 
called grid A) with masses in the range 0.42 - 0.47 M in steps 
of 0.01 M . The maximum mass we obtained for the degenerate 
core of an RGB star is 0.47 M , thus we did not consider post- 
flash sdB stars above this mass. The H-envelope masses, M env ,o, 
considered are 0.0001, 0.0003, and 0.0006 M , where M enVi0 is 
defined as the total mass of the hydrogen content directly af- 
ter the removal of the envelope. Thus, we have 18 sdB evolu- 
tion tracks, which we followed until the end of He-core burning. 
After each 10 7 yr of sdB evolution, the seismic properties were 
calculated, and we only considered models with unstable modes. 
We have not found unstable modes in our post-He-core burning 
models. Therefore our analysis is limited to sdB stars in their 
He-core burning phase, leading to a total of 402 seismic models 
in the range of T eS = 25, 000 - 34, 000 K and log g = 5.4 - 6.0. 
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Fig. 2. Evolutionary tracks in the log g— T e ff diagram of sdB stars 
with different iron profiles. The three tracks start with the same 
zero-age EHB model with M sdB = 0.47 M G and M env ,o = 10~ 4 
Mo, created from a 1.00 M Q ZAMS model. The lowest solid 
curve is for a model with no iron enhancement, the middle dot- 
ted curve is for a Gaussian iron increase centered at log T = 5.3, 
and the upper dashed curve is for iron increased with a factor 10 
uniformly in the envelope. The labels I, II, II indicate the ZAHB, 
the end of core-He burning and the end of He-shell burning, re- 
spectively. It takes the sdB star 1.8 x 10 8 yr to evolve from I to 
II, and 10 7 yr from II to III. 




Fig. 3. Left: the frequencies of the stable (x) and unstable (+) 
modes (Z < 2) as a function of the EHB age, for the three tracks 
given in Fig. [2] Right: the opacity as a function of the temper- 
ature after 10 yr of sdB evolution. Panels (a&d) correspond to 
/ = 10, (b&e) are for the Gaussian parametrization of /, and 
(c&f) are for the model with no iron enhancement. 



The grid of non-canonical sdB stars (hereafter called grid B) 
consists of 5 tracks: (M sdB (M Q ), M env o (M )) = (0.44, 0.005), 
(0.45, 0.005), (0.46, 0.005 ), (0.47, 0.0075 ), and (0.47, 0.005). 
Along these tracks we have in total 98 seismic models, again 
taken after each 10 7 yr of sdB evolution. In Table[T] more details 
about the models are given. 

The sdB evolution tracks and the seismic models can be seen 
in Fig. [4] The tracks start directly after the removal of the enve- 
lope. For the post-flash models, this corresponds to the zero-age 
EHB. The post-non-degenerate models have hydrogen extending 
to deeper layers, hence allowing some H-shell-burning, before 
reaching the zero-age EHB. Note that, although the two differ- 
ent types of sdB stars can evolve through the same logg - T e s 
during core-He-burning, the post-He-core-burning evolution dif- 
fers, see Fig. The post-non-degenerate sdB star again has a 
short phase (~10 6 yr) of H-shell burning, before starting He- 
shell burning. 

3.3. Comparing two representative models 

We examined the physical differences in the interior structure of 
a post-flash (a) and a post-non-degenerate (B) sdB with same 
logg and T e ff. We chose as representative the models circled in 
Fig. [4}} at logg = 5.78 and T e s = 30 kK. One of the main differ- 
ences is the abundance profiles, see Fig. [6^. The He-H transition 
layer of B is much broader and located deeper in the star. The 
envelope of B is in the region where the shrinking convective 
core passed through during the MS, hence the low H-abundance 
here: X — 0.18. For a, the helium core has grown into the region 
that used to be part of the convective envelope during the RGB, 
and, as a result, the He-H transition region is much narrower 



while the H-abundance in the envelope is around X = 0.66. We 
will discuss the possibility of diffusion of hydrogen to the sur- 
face in ^4] In Fig. I5H, we show the iron mass fractions, and in 
Fig. [6^ the resulting opacity profiles. The outer opacity bump 
near log T = 4.7 is associated with helium ionization, and the 
iron opacity bump near log T = 5.3 is enhanced by the local iron 
accumulation. 

Two important quantities in stellar pulsation theory are the 
the Brunt-Vaisala frequency N and the Lamb frequencies L/, 



N z 



Gm 6 

~J2 



ad ' 



v + %v„ 



(4) 



<51np dlnp 
with o = — — — — and <p = 



d\nT 



din// 



1(1 + \)ct 



(5) 



where n is the molecular weight and c s is the adiabatic sound 
speed. When N 2 < 0, the Ledoux criterion for dynamical stabil- 
ity is violated fLedoux 1947[ ). Thus, in Fig.|6]5&c, the convective 
regions can be clearly identified. The innermost one is related to 
the convective core, and the outer two are the narrow convective 
layers due to iron (near log T = 5.3) and helium ionization (near 
log T = 4.6). Also, chemical gradients are apparent in A^ 2 in the 
form of localized peaks. The peak near the center is identified 
with the C-O/He transition zone and the outermost peak to the 
He-H transition zone. The Lamb frequencies are plotted in the 
same figure to indicate the propagation zones of the g-modes, 



Hu et al.: A seismic approach to testing different formation channels of subdwarf B stars 



5 






5.4 




O.O 




5.6 


CM 


5.7 


CO 




E 






5.8 








5.9 


O 




6 




6.1 




6.2 



24 



(b) 



H-shell 



H-shell 





65 



60 



55 



50 45 
I,,,, (1000 K) 



40 



35 



30 



Fig. 4. (a) sdB Evolutionary tracks in the log g - T e ff diagram. The solid curves correspond to the post-flash sdB models, and the 
dotted curves to the post-non-degenerate sdB models. The circles and squares indicate models on which we performed pulsation 
calculations for the post-flash and post-non-degenerate tracks respectively, (b) Evolutionary tracks of a post-flash sdB star created 
from of low-mass (Mzams = 1 -00 M ), and a post-non-degenerate sdB star created from a high mass (Mzams = 3.00 M Q ) progenitor, 
given by the solid and dotted curve respectively. Along the tracks the dominant energy source (either H-shell-burning, He-core- 
buming or He-shell-burning) is noted. The bold part of the tracks indicate the He-core-burning phase and the thin part the He-shell 
burning phase. During the dot-dashed parts of the post-non-degenerate track, H-shell burning is the dominant energy source. At the 
circle in (b), we selected from both tracks a model for a detailed comparison, see Section 3.3 



Table 1. The models in grid A and B.^\ 



track 


# models 


M sdB (Mq) 


M ell v,0 (Mq) 


M Z A 


Al 


27 


0.42 


0.0001 


1.85 


A2 


27 


0.42 


0.0003 


1.85 


A3 


27 


0.42 


0.0006 


1.85 


A4 


25 


0.43 


0.0001 


1.80 


A5 


25 


0.43 


0.0003 


1.80 


A6 


25 


0.43 


0.0006 


1.80 


Al 


23 


0.44 


0.0001 


1.75 


AS 


23 


0.44 


0.0003 


1.75 


A9 


23 


0.44 


0.0006 


1.75 


A10 


21 


0.45 


0.0001 


1.65 


All 


21 


0.45 


0.0003 


1.65 


A12 


21 


0.45 


0.0006 


1.65 


A13 


20 


0.46 


0.0001 


1.55 


A14 


20 


0.46 


0.0003 


1.55 


A15 


20 


0.46 


0.0006 


1.55 


A16 


18 


0.47 


0.0001 


1.00 


A17 


18 


0.47 


0.0003 


1.00 


A18 


18 


0.47 


0.0006 


1.00 


SI 


23 


0.44 


0.005 


2.75 


52 


18 


0.45 


0.0075 


2.75 


53 


21 


0.46 


0.005 


2.90 


64 


18 


0.47 


0.0075 


2.90 


B5 


18 


0.47 


0.005 


3.00 



" The label of the track given in the first column is used for reference 
in Fig. [7] The second column gives the number of seismic models along 
that track. The final column shows the mass of the ZAMS models from 
which the sdB star is created. 



cr 2 < (N 2 ,L 2 ), and the p-modes, cr 2 > (N 2 ,L 2 ), where cr is the 
angular pulsation frequency. It is apparent that g-modes are deep 
interior modes, while p-modes probe the superficial outer layers 
as pointed out by Charpin et et al.| (pOOO). 



Clearly, models a and f3 have very different physical char- 
acteristics. To establish how this affects their seismic properties, 
we compare their frequencies in Fig. [5] Since the large frequency 
separation A/ = f„j - f n -\,i is mainly dependent on the dynam- 
ical timescale, we see that Af at high frequencies is more or 
less the same for the two models. The lower frequencies, how- 
ever, are in fact mixed modes that are more sensitive to the core, 
and we see a better distinction between models a and ft here. 
Moreover, the frequency ranges of excited I — - 2 modes 
are not the same for these two models; for model a this range 
is [3.4 mHz, 10.6 mHz], and for model /?, it is [5.7 mHz,14.2 
mHz] . The excited modes thus have lower frequencies in model 
a than in model ft. To understand this, we compare in Fig.[6ji the 
work integral of these two models for the radial mode pj. The 
work integral increases towards the surface in the driving region 
and decreases towards the surface in damping regions. The sur- 
face value is the dimensionless growth rate, positive for unstable 
modes and negative for stable ones. We see that p-j is unstable in 
model ft, but stable in model a, which can also be seen in Fig. [5] 

A first possible origin of the differences could come from the 
opacity, since the driving is a K-mechanism operating in the iron 
opacity bump. Fig. [6^ shows that the opacity is slightly larger 
for model ft. The driving is thus a little more efficient in model 
ft. But this is not the main source of differences. Since the enve- 
lope H-fraction is much smaller in model ft (X — 0.18) than in 
model a (X - 0.66), the molecular weight is larger and, at given 
temperature, the density is significantly higher (~1.5x) in model 
ft, as shown by Fig. |6f. The driving of the modes is related to 
the opacity, which is mainly a function of temperature. Hence, 
if the eigenfunctions of two given modes have the same shape 
as a function of temperature, the driving is the same. Here we 
compare the modes pi of two models with the same radius. As 
it is usually found for p-modes, their last node is located at the 
same geometrical distance from the surface Ar. But the gradient 
of temperature is not the same for the two models: |dT/dr| oc Kp 
is greater in model ft because of the higher density, as can be 
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Table 2. The models with minimum Ai 2 for the four different 
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where // is one of the rib excited frequencies of star b, and f' a 
is the correspondingly matched frequency of star a, expressed 
in mHz. The frequency matching is done such that Ai 2 
imized by brute-force fitting. It is clear that the lower Ai 



Fig. 5. The pulsation frequencies for modes / = - 2. Solid line 
segments are for model a and dotted ones for model /3. Unstable 
modes are given by long line segments, and stable modes by 
short line segments. The modes with radial order n — 1 and 
radial mode p 7 are indicated. 



seen in Fig^&g. Hence, the difference between the tempera- 
ture at the last node and at the surface, \AT\ |d7*/dr| Ar, is 
greater in model B than in model a. This is exactly what we find 
in Fig. [6}, where the eigenfunction \6T/T\ is given. In terms of 
the temperature, the last node is closer to the surface in model a 
than it is in model B. To get the same driving as in model B, the 
last node of model a would have to be deeper in the star, which is 
only possible by considering a mode with lower radial order and 
frequency. Hence, the frequencies of excited modes are lower in 
model a than in model B. 

3.4. Comparing two grids of models 

We investigated if it is possible to distinguish between a post-He- 
flash (a) and a post-non-degenerate sdB stellar model (b) from 
observed oscillation modes. Imagine we observed the frequen- 
cies of b; is it then possible to find an acceptable seismic match 
in our grid of canonical post-flash models (grid A)? We took as 
'observed' frequencies those of unstable modes up to / = 2. We 
did this for each model b in grid B, thus finding the best seismic 
match within grids A and B. 

Since frequency separations follow from asymptotic rela- 
tions for p-modes, the frequency is a natural quantity for model 
comparison. Despite this, periods have been used more often in 
the literature so far, when comparing observed modes of sdB 
stars with those predicted by models. We also considered pe- 
riod matching, but found frequency matching more suitable to 
compare the p-modes of the models. This will be different for 
g-modes, where the mode period is the natural quantity to com- 
paring observations with models. 

To quantify 'acceptable', we used the merit function 
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better the match between a and b. Matches with M 2 > 0.05 are 
considered unacceptable, which is a generous limit, as we will 
see later. We investigated four different scenarios: 

(i) We are not able to identify the modes, log g and T e ff of the 
'observed' star b are unknown, and the 'observed' frequen- 
cies are allowed to be matched with both stable and unstable 
frequencies of the 'theoretical' model a. 

(ii) Same as (i), except the modes are identified, thus the /-value 
must be matched. 

(iii) Same as (i), except logg and T & g are known within errors of 
dlogg = 0.1 and dT eff = 1000Z. 

(iv) Same as (i), except the 'observed' frequencies are only 
matched with unstable 'theoretical' frequencies, i.e. assum- 
ing that the theory correctly predicts which frequencies are 
exited and which are not. 

In Fig. [7] we show M 2 for each gridpoint in grids A and B 
for the scenarios (i)-(iv). The matches with low Ai 2 are visible as 
dark diagonal regions, This is a result of the change in frequen- 
cies during the sdB evolution. From Fig.[7Ji)-(iii), it is clear that 
the distinction between models a and b is drastically increased 
if we have either mode identification or spectroscopic logg and 
T e ff values. Fig. |7jiv) shows that, if we only allow matching to 
unstable (and not stable) frequencies of a, there are no matches. 

The matches with lowest Ai are circled in Fig. |7]and details 
of these models are shown in Fig.|8]and Table|2] For all scenarios 
the same model b gives the best match, namely the last model 
of evolutionary track B5. We understand that in terms of only 
the higher frequencies with radial order n > 3 being excited in 
this model. As we discussed in j |3.3| the lower frequencies are 
more sensitive to the deeper layers, thus the distinction between 
models a and b is better detected at low frequencies. 

As a comparison to a real case, we considered the optimal 
model for PG 0014+067, for which |Brassard et al-lpOOT) found 
X 2 = 0.5374, where x 2 is a merit function based on mode period 
comparison. Translated to our frequency merit function, this is 
equivalent to Ai 2 = 0.0084. Although we find, in principle, seis- 
mic matches between a and b with Ai 2 of this order for scenarios 
(i)-(iii), they are not statistically favoured. For scenario (i) we 
find that 12 of the 98 models in grid B can be matched with a 
model in grid A with Ai 2 < 0.01, and this is 7 for scenario (ii), 
and only 1 for scenario (iii). 
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Fig. 6. Physical and seismic quantities of the two representative models a (post-flash) and/? (post-non-degenerate). The profiles are 
shown as a function of the temperature: (a) the hydrogen (X) and helium (Y) mass fractions; (b)&(c) the Brunt- Vaisala and Lamb 
frequencies for model a and /3, respectively; (d) the iron mass fraction; (e) the opacity; (f) the density, we also plotted the density 
ratio pplp a on the right axis; (g) the temperature gradient \dT/dr\; (h) the work integral for radial mode pi; and (i) the eigenfunction 
\5T/T\ for the radial mode p-j. In all panels except (b)&(c), solid lines represent model a and dotted lines model /3. Note that panels 
(a)-(c) give the profiles throughout the entire star, while the profiles in panels (d)-(h) are for the stellar envelope. 



4. Discussion & conclusions 

We studied the so far neglected, post-non-degenerate sdB stars 
and compared their physical and seismic characteristics with 
those of canonical post-flash sdB stars, both formed in the CEE 
channel. The results presented here are a first step in distinguish- 
ing these two kinds of sdB stars on the basis of their observed os- 
cillation character, which is necessary if seismic modelling is to 
achieve reliable mass determination. Furthermore, the observa- 
tion of a post-non-degenerate sdB star in a post-CE binary would 
give strong constraints on the CE evolution. We plan to continue 
such investigations with an application to the sdB pulsator in the 
post-CE, eclipsing binary PG 1336-018 which we started in Hu 
[etaT1 ( [2D07l > and | Vuckovic etaL] ( |2007) . 

We find that, in principle, a post-non-degenerate sdB star 
may appear as an EC 14026 star with similar pulsation frequen- 
cies as the canonical post-He-flash sdB star, although it is not 
likely. Additional observables, such as spectroscopic \ogg and 
r e ff determinations and/or empirical mode identification from 
observables enable us to distinguish the two types of sdB stars 
more decisively. The frequency range of the unstable modes is 
also an important discriminator between the two formation chan- 
nels. In general, for the same logg and r e ff values, the excited 



frequencies of the post-non-degenerate sdB star are higher than 
the excited frequencies of the post-flash star. This is a direct re- 
sult of the differing interior structures. Thus, special attention 
must be paid when observed frequencies are matched with theo- 
retically predicted ones of modes that are not excited. 

Up to now, there have not been any evolutionary models of 
sdB stars available that include the coupling between diffusion 
and evolution consistently. This is a deficiency, since iron accu- 
mulation due to radiative levitation is responsible for the pul- 
sational instability in these stars ( Charpinet et al.jfl 996| >. Also, 
it has been shown by Fontaine et al. ( 2006} thalthe iron accu- 
mulation changes the frequencies significantly. In our study, we 
have parametrized the iron accumulation, so that we can, at least 
in an approximative manner, simultaneously take into account 
the effects of iron enhancement and evolution on the pulsation 
modes. 

Here we have not considered the influence of the other dif- 
fusive processes, i.e. diffusion due to gradients of pressure, tem- 
perature, and concentration. To a certain extent this can affect 
our results, because one of the main differences between the 
two types of sdB stars is the chemical composition of the stel- 
lar envelope. Specifically, we find in the envelope of the post- 
non-degenerate sdB star an H-mass fraction of X = 0.18, while 
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Fig. 7. Colourmap of Ai 2 . The models of grid A are plotted along the vertical axis, and the models of grid B along the horizontal 
axis. The models are divided in blocks according to evolutionary track, where we have ordered the blocks with increasing mass 
(see Table|2|, and in the blocks the models are ordered with increasing age. The four panels are for the four different scenarios (i), 
(ii), (ii), and (iv), as described in the text. In each panel with have pointed out and circled the gridpoint with minimum Ai 2 . The 
frequency matching of these gridpoints can be seen in Fig. [8] Note that panel (iv) has no acceptable matches, but we still circled the 
one with the lowest M 2 . 
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Fig. 8. The frequency matches of the minimum Ai 2 gridpoints. 
The /-value of model a is depicted on the vertical axis. For sce- 
narios (i) and (ii), this is equal to the /-value of model b, for 
scenario (iv) it is not. Solid line segments are for model a and 
dotted ones for model b. In the upper panel, we have indicated 
the lowest radial order of the unstable modes of model b. 



the post-flash sdB star has X = 0.66 there. Normally, it is as- 
sumed that sdB stars have H-rich or even pure H-envelopes, 
caused by gravitational settling. While this is true for the out- 
ermost layers, diffusion is not expected to work efficiently at 
depths log 4 > -3 ( |Richard et aLp002| pichaud et al.||2057) . 
Since the envelopes of the post-non-degenerate sdB stars ex- 
tend to logq > -2 (i.e. T > 10 7 K), we do not expect diffu- 
sion to wash away all the qualitative differences in the chemi- 
cal profiles, although the differences may be less pronounced. 
Diffusion, however, will significantly change the surface abun- 
dances of our models and likely will bring them into agreement 
with the observed values. Spectroscopic line profile analysis has 
shown that the majority of sdB is He-deficient, and only a few 
are He-rich ( jEdelmann eTaLl|2003l |Lisker et al.l[2004] >. Stellar 



evolution models that include diffusion coupled to reliable at- 
mosphere models are needed to assess whether the two different 
formation channels will be distinguishable via a spectroscopic 
abundance analysis. We are currently computing such evolution- 
ary sdB models including diffusion due to gradients of pressure, 
temperature, and concentration. Our preliminary results indeed 
agree with our expectations, i.e. the H-surface abundance in- 
creases on a very short timescale, but the chemical profiles at 
deeper layers are not affected. The pulsational properties of these 
improved models will be discussed in detail in a forthcoming pa- 
per. 

We have made a modest grid of models that is sufficient 
for our comparative study. Detailed seismic modelling of an ob- 
served star, however, will require a finer grid. For now, we have 
chosen not to make sdB models above 0.47 M Q , since this is 
the maximum mass the degenerate He-core of a red giant with 
Z = 0.02 can have before experiencing the He-flash. A metal- 
licity of Z = 0.004 allows the He-core to grow up to 0.48 M Q 
on the RGB. However, we find that, in order to excite modes in 
these low metallicity stars, an iron enhancement greater than a 



factor 10 is required. This was to be expected, since Charpinet 
et al. ( |1996 i found unstable pulsation modes for models with 
uniform Z > 0.04 in the H-rich envelope. We have therefore not 
pursued these models further. The question whether post-flash 
sdB stars can have masses > 0.47 M Q is also closely related to 
the input physics (e.g. convective overshooting) and the physics 
of the He-flash, and needs to be examined further. 

In this paper, we have focused on the short-period p-mode 
sdB pulsators. The case of the long-period g-mode sdB pulsators 
is, although challenging from an observational point of view, an 
additional very interesting theoretical case study. The p-modes 
only probe the outermost layers, and hence are less affected by 
the differing composition gradients than the g-modes, as they 
propagate deeper into the star. The long-period sdB pulsators are 
interpreted as cooler sdB models with much thicker hydrogen- 
envelopes than the short-period sdB pulsators (Fontaine et al.| 
2003[|Jeffery & Saio 2006 ). Since the g-modes are deep interior 
modes, full evolutionary models including iron accumulation, as 
developed here, are required to model these stars. At present, 
these are not available yet. We are currently developing a similar 
approach to the one presented here to study the long-period sdB 
pulsators. 
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